**# Bookmark #1 Setting

*Windows
global root  = "G:/Dropbox/Environmental Injustice/Empirical"
global rawdata = "$root/rawdata"
global workdata = "$root/workdata"
global results = "$root/results"
global outfile = "$root/outfile"

cd "$workdata"


**# Bookmark #2 Regression
*Cardiovascular Disease
use "$outfile/Table S21 heartandday.dta",clear

reghdfe heart lowedu5_post_notapwater lowedu5_post lowedu5_notapwater notapwater_post lowedu5 notapwater post pm25 so2,absorb(i.id i.province#i.year) cluster(i.province i.year)

outreg2 using "$results/Table S21.xls",replace ctitle(Cardiovascular Disease) ///
keep(lowedu5_post_notapwater lowedu5_post) ///
addtext(Individual FE, YES, ProvincebyYear FE, YES) ///
bdec(3) tdec(3) rdec(3)

*Inactive Days
reghdfe day lowedu5_post_notapwater lowedu5_post lowedu5_notapwater notapwater_post lowedu5 notapwater post pm25 so2,absorb(i.id i.province#i.year) cluster(i.province i.year)

outreg2 using "$results/Table S21.xls",append ctitle(Inactive Days) ///
keep(lowedu5_post_notapwater lowedu5_post) ///
addtext(Individual FE, YES, ProvincebyYear FE, YES) ///
bdec(3) tdec(3) rdec(3)

*Tumor
use "$outfile/Table S21 tumor.dta",clear

reghdfe tumor lowedu5_post_notapwater lowedu5_post lowedu5_notapwater notapwater_post lowedu5 notapwater post  pm25 so2,absorb(i.id i.province#i.year) cluster(i.province i.year)

outreg2 using "$results/Table S21.xls",append ctitle(Tumor) ///
keep(lowedu5_post_notapwater lowedu5_post) ///
addtext(Individual FE, YES, ProvincebyYear FE, YES) ///
bdec(3) tdec(3) rdec(3)

*Daily Work Hours
use "$outfile/Table S21 hour.dta",clear

reghdfe hour lowedu5_post_notapwater lowedu5_post lowedu5_notapwater notapwater_post lowedu5 notapwater post  pm25 so2,absorb(i.id i.province#i.year) cluster(i.province i.year)
outreg2 using "$results/Table S21.xls",append ctitle(Daily Work Hours) ///
keep(lowedu5_post_notapwater lowedu5_post) ///
addtext(Individual FE, YES, ProvincebyYear FE, YES) ///
bdec(3) tdec(3) rdec(3)

*Wage_100
use "$outfile/Table S21 wage.dta",clear

reghdfe wage_100 lowedu5_post_notapwater lowedu5_post lowedu5_notapwater notapwater_post lowedu5 notapwater post  pm25 so2,absorb(i.id i.province#i.year) cluster(i.province i.year)
outreg2 using "$results/Table S21.xls",append ctitle(Wage_100) ///
keep(lowedu5_post_notapwater lowedu5_post) ///
addtext(Individual FE, YES, ProvincebyYear FE, YES) ///
bdec(3) tdec(3) rdec(3)
	
	






